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This talk summarizes the results of a phenomenological analysis of heavy quark fragmentation data published by the 
CLEO and BELLE collaborations at yfs = 10.6 GeV and by the LEP collaborations at i/s = 91.2 GeV. Several 
theoretical ingredients are employed: next-to-leading order initial conditions, evolution and coefficient functions; soft- 
gluon resummation to next-to-leading-log accuracy; a next-to-leading order matching condition for the crossing of 
the bottom threshold in the evolution. Important initial-state electromagnetic radiation effects in the CLEO and 
BELLE data are also accounted for. We find that with reasonably simple choices of a non-perturbative correction to 
the fixed-order initial condition for the evolution, the data from CLEO and BELLE can be fitted with remarkable 
accuracy. The fitted fragmentation function, when evolved to LEP energies, does not however represent fairly the D* 
fragmentation spectrum measured by ALEPH. Large non-perturbative corrections to the coefficient functions of the 
meson spectrum are needed in order to reconcile CLEO/BELLE and ALEPH results. 



1. Introduction 

The CLEO [1] and the BELLE [2] Collaborations have recently published high-statistics and high-accuracy data 
for various charmed mesons fragmentation in e + e _ collisions at a centre-of-mass energy of 10.6 GeV. These data are 
the first ones able to rival in quality (and, in fact, to best) similar ones published a few years ago by the ALEPH 
Collaboration [3] at 91.2 GeV. 

Taken together these sets of data allow for a phenomenological analysis [4] which spans a fairly large energy gap. 
It is therefore possible not only to test the ability of the theoretical framework to describe the data well, but also to 
perform the evolution from one energy to the other, and look for evidence (or absence) of power suppressed effects. 

The theoretical framework employed is based on a next-to-leading order QCD description of the fragmentation 
process [5] , including a next-to-leading log accurate resummation of collinear and soft gluons [6] 1 . 

A few more details complete the theoretical picture. They are briefly summarized below, and they are outlined in 
detail in [4]. 

First, and most important, the soft-gluon resummation needs to be artificially regularized at large x (or, equiva- 
lently, large N in moment space). The choice of the regularization defines the perturbative distribution and therefore 
directly influences the non-perturbative one extracted by fitting the experimental data. In this work we have decided 
to use a prescription which, while being as simple as possible, possesses the following desirable features: 

(i) it is consistent with all known perturbative results, 

(ii) it yields physically acceptable results, 

(iii) it does not introduce power corrections larger than generally expected for the processes in question, i.e. NA/m 
for the initial condition [9-12] and NA 2 /q 2 for the coefficient functions [13], where A is a typical hadronic scale 
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1 It is worth noting that very recently next-to-next-to leading order results for the time-like non-singlet splitting function have become 
available [7]. Together with the previously available Q[a%) initial conditions [8] they make possible to repeat at least part of this analysis 
with higher accuracy. 
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of a few hundreds MeV. 

In practical terms, our choice will yield a fragmentation function which does not become negative in the large- A 
region. This will allow for a good description of the data up to the the x = 1 cndpoint. 

Second, while evolving the charm fragmentation function through the bottom threshold, one needs in principle to 
modify the evolution equations and properly match at this threshold. Moreover, production of charm via gluon split- 
ting should be allowed, since it represents a non entirely negligible source at LEP energies. Both these features have 
been implemented in [4], thus departing from the simple non-singlet only description previously usually employed. 

Finally, the experimental data as measured by CLEO and BELLE still contain the effect of electromagnetic initial 
state radiation. We have estimated that this effect is not negligible in this case (as it is, instead, at LEP, due to the 
physical cutoff provided by the Z Q resonance peak) . We have therefore proceeded to simulate it and to deconvolute 
it from the data before fitting them with a pure QCD description of the fragmentation process. Again, details can 
be found in [4]. 



2. Non-perturbative fragmentation function 

In the heavy-quark fragmentation-function formalism, the largest non-perturbative effects come from the initial 
condition, since one expects power corrections of the form A/m. We assume that all these effects can be described by 
a non-perturbative fragmentation function Z?np, that takes into account all low-energy effects, including the process 
of the heavy quark turning into a heavy-flavoured hadron, that has to be convoluted with the perturbative cross 
section. Thus, the Mellin transform of the full resummed cross section, including non-perturbative corrections, is 



a H (N, q 2 ) = a Q (N, q 2 ,m 2 )D NP (N) 



(1) 



We have attempted to fit CLEO and BELLE D* data using several forms for Z?np- We found that the best fits are 
obtained with the two-component form 



I 



D NP {x) = Norm, x S(l - x) + cN~l(l - x) a x l 



(2) 



with 



N a , b 



Jo 



1 



xrx" 



(3) 



This form is a superposition of a maximally hard component (i.e. the delta function) and the form proposed in 
Ref. [14]. It can be given a simple phenomenological interpretation, the hard term corresponding in some sense to 
the direct exclusive production of the D*, and the Colangelo-Nason form accounting for D*'s produced in the decay 
chain of higher resonances. 

Following the approach of Ref. [15], we assume that the D meson non-perturbative fragmentation function is the 
sum of a direct component, which is isospin invariant, plus the component arising from the D* decay. The decay 
D* —> Dir is very close to threshold, so that the D has the same velocity of the D*, and their momenta are thus 
proportional to their masses. Under these circumstances, the component of the D fragmentation function arising 
from D* — > Dir decays is given by 



B(D* Dir) D°{x) 



(4) 



where we have defined 



D* 
NP 



n%D 



m D 



mp 



(5) 
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Eq. (2): a = 1.8 ± 0.2, b = 11.3 ± 0.6, c = 2.46 ± 0.07, total x 2 = 139 


Set 


(C) D*+ 


(B) D*+ -» D° 


(B) ZT+ -> D+ 


(C) D*° 


(B) D*° 


Norm. 


0.238 


0.253 


0.227 


0.225 


0.211 


X 2 Ms 


33/16 


63/46 


13/46 


13/16 


17/46 



Table I: Results of the fit to D* CLEO (C) and BELLE (B) data. The last line reports the \ 2 over the number of fitted points 
for each data set. 



and B(D* — ► Dn) is the branching ratio of D* — ► £>7r. Observe that has been defined so as to have the same 
normalization as D® P . In N space we obtain immediately 



dZ(n) = d°;(n) 



m D 



-i N-l 



(6) 



For the D* — ► D7 decay, in the D* frame, the D has non-negligible velocity, but it is non-relativistic, its momentum 
being given by 



Pd 



It can easily be shown [4] that, in moment space, we can write 

m D * (m D +Pd) N - (m D - Pd) N 



D°(N)=D°;(N) 



2 P d 



(7) 



(8) 



We thus describe D + /° production as the sum of a primary (i.e. not coming from D* decays) component, plus the 
contributions coming from D* decays 



D&{x) = D°;*{ X ) + B{D*+^D+^)b^{x) 

+B(D*+ -» D+-y)D° + D(x) , 
D^(x) = D°° p - P (x) + [B(D*+ -> D a 7r+) + B(D*° -> D\°)}D°° (x) 
+B(D*° — > D -f)D®° (x) . 



(9) 
(10) 



3. D mesons data fits near the T(45) 

Several parameters enter our calculations. First of all, at all matching points, there are scale choices that could be 
varied, to yield a perturbative uncertainty in our result. Those are the initial evolution scale /zo, the matching scale 
for the crossing of the b threshold /i t hr, and the final evolution scale \x. In the present work we fix 

H =m, [1= \[q* = yfs , jlithr = TOthr = m b . (11) 

These scales could, in principle, be varied by a factor of order two around the values listed above, yielding a sensibly 
different result. However, in general, the scale variation will simply result in different values for the fitted parameters 
of the non-perturbative form. When computing cross sections for different processes, one should then use the 
parametrization appropriate for the scale choice that has been made in the fit, hence compensating for the change. 
In the present work we will not pursue this issue further, since our aim is simply to show that a fit within QCD is 
possible. A similar remark applies to the value of A QCD and the quark masses, that we will fix at 

Aqq D = 0.226 GeV , m c = 1.5 GcV , m b = 4.75 GeV . (12) 

The result of the fit is reported in Table I, and in Figs. 1, 2 we show some of the data and the corresponding fitted 
curve, both in x and moment space. 
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Figure 2: Fit to BELLE D* + -> D° data. 



4. D mesons data fits on the Z° 

Fig. 3 shows a similar fit to ALEPH D* + data [3]. We fit the data in the region x £ [0.4, 1] using the non-singlet 
component only, since a subtraction of the gluon-splitting contributions was performed by ALEPH. Observe that, in 
this calculation, the bottom-threshold crossing has to be dealt with. We also show, for comparison, the full evolution 
result (dashed line), using the same parameters obtained in the non-singlet fit. As we can see, the difference is only 
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visible at small x. The result of the fit for the non-perturbative parameters is 

a = 2.4±1.2, 6 = 13.9 ± 5.7 c = 5.9±1.7, (13) 

with a x 2 — 4-2 for 13 fitted points. These results are not really consistent with those for the T(45) data in Tab. I. 

In order to better quantify the discrepancy between Eq. (13) and Tab. I we use the parametrization of CLEO and 
BELLE data to predict the D* fragmentation function at LEP energies. The LEP prediction, using the parametriza- 
tion of Table I, is reported in Fig. 4 together with the ALEPH data. We find a x 2 = 60.1 (for 13 fitted points) for 
this parametrization. Thus, the description is not satisfactory, especially in the large- x (large- N) region. 

In Fig. 5 we show the ratio of the moments of ALEPH D* + data over our prediction. We observe that the N 
dependence of the ratio is well described by the functional form 

(14) 



1 + 0.044 (N — 1) 



where, since the first moment of the non-singlet distribution should be exactly given by the theory (because of charge 
conservation), we normalize to one the extrapolation of the data to N = 1. 

We can only speculate about the possible origin of the discrepancy and the form of the coefficient of (N — 1) in 
Eq. (14). Assuming that we are dealing with a non-perturbative correction to the coefficient function of the form 

1 + <») 



this would lead to the extra factor 

1 4- C(N-l) 

(16) 



1 , C(N-l) 

1 + Mt 



(where Mz and My are the Z° and T(45) mass) to be applied to our prediction for the ALEPH data. For C = 5 GeV 2 
we reproduce the behaviour of Eq. (14). In Ref. [13], on the basis of a calculation of infra-red renormalon effects, 
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Figure 4: ALEPH D* + data, compared to the QCD prediction. 



a 1/q 2 power correction is found, with an TV dependence marginally compatible with (15). No 1/yY correction is 
found. Ref. [16] also predicts a leading 1/q 2 power correction. However, the C w 5 GeV 2 coefficient would appear to 
be somewhat too large 2 . Alternatively, if we admitted the existence of corrections to the coefficient functions of the 
form 



1 , C(N-l) 

then we would find C s=s 0.52 GeV, a much more acceptable value. We observe that a form 

N-2 



(17) 



1 + 



C 



1 + 



C(N-2) 



was required in Ref. [17] to fit light-hadron fragmentation data. 



(18) 



Demonstrating the absence (or the existence) of l/\/q 2 corrections in fragmentation functions would be a very 
interesting result, since it would help to validate or disprove renormalon-based predictions. Unfortunately, the low 
precision of the available data does not allow, for the time being, to resolve this issue. 

We would like to remark that the discrepancy between the CLEO/BELLE and ALEPH data exclusively depends 
upon the evolution between the T(45) and Z° energies. The method we used to describe the CLEO/BELLE data 
(i.e. the perturbative calculation of the fragmentation function, the Sudakov effects in the initial conditions and the 
parametrization of the non-perturbative part) does not affect the conclusions of the present section. In fact, we 
can simply compute the ratio of the moments of the inclusive D* + (ISR corrected) distribution at CLEO/BELLE 
and ALEPH, and compare it to the theoretical prediction. The result of this comparison (where we have used, for 
simplicity, BELLE data only) is displayed in Fig. 6. The curves are given by 



2 If we believe that it is the maximum meson energy, not y/lp, that controls power effects, than we would have C ~ 1 GeV 2 , a more 
acceptable value. 
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Figure 5: ALEPH D* + data, compared to the QCD prediction. 
a Q (N,M 2 ,m 2 ) _ a q (N,M 2 . /f i 2 z ) 2 , l + a s (^)/n 

a Q (N,M 2 ,m 2 ) l + a s ( M |)A 1 ' ^' Mx ' a q (N, M 2 , [U) 

where [iz and fix are the factorization scales and the evolution factor E is given by the solution of the Altarelli-Parisi 
evolution equation. Notice that low-scale effects, both at the heavy quark mass scale and at the non-perturbative 
level, cancel completely in this ratio, making its prediction entirely perturbative. As we can see from the figure, 
the rather large scale uncertainty displayed by the NLO result is much reduced when Sudakov effects are included. 
In both cases, however, the data clearly undershoot the pure QCD prediction, being instead compatible with the 
inclusion of the correction factor (14) (dotted lines). 

One can legitimately wonder whether some of the theoretical ingredients might hide further uncertainties. We have 
checked that the regularization procedure needed to deal with the Landau pole in the soft-gluon resummation has 
very little impact on our curves. Using the very large value Aq 5 c D = 0.3 GeV would lower the theoretical predictions 
by no more than 11% for N < 20, very far from explaining the observed effect. The deconvolution of ISR effects, 
that hardens the T(AS) data, but is insignificant on the Z , widens the discrepancy, which would however still be 
partially visible even if the data were not corrected for e.m. radiation. 

Because of the relatively low energy of the data on the T(4S I ), it is also legitimate to wonder whether charm-mass 
effects could be responsible for the discrepancy between LEP and T(45*) data. We have not included mass effects in 
the present calculation. However, in Ref. [18], mass effects in charm production on the T(45) where computed at 
order a 2 , and found to be small. In fact, they amount to a correction of the order of 1% at N — 5, 4% at N = 11 
and 7% at N = 20, very far from being able to explain our observation. We thus believe that it is unlikely that mass 
effects could play an important role in explaining this discrepancy 3 . 



3 In Refs. [10] and [9], on the basis of the analogy with the spacelike case, corrections of the form Am/9 2 are introduced. The importance 
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N 

Figure 6: The ratio of ALEPH and BELLE moments for the D* + fragmentation function, compared to QCD evolution. The 
solid band is obtained with QCD NLO evolution and Sudakov effects in the coefficient functions, while the dashed bands is 
NLO evolution only. The bands are obtained by setting fi z /r = £M z /r and varying 1/2 < £ < 2. 

5. Conclusions 

This phenomenological analysis [4] of heavy quark fragmentation in e + e~ collisions shows that it is possible to 
perform excellent fits of D* /D meson fragmentation spectra in perturbative QCD, using all known results on the 
perturbative heavy-quark fragmentation function, and compounding them with a simple parametrization of non- 
perturbative effects. 

A second striking result is the evidence of large non-perturbative effects, visible in the relation between the 
D* fragmentation function at the T(AS) and Z° energies. It would be interesting to understand the power law 
of these contributions. Their magnitude would suggest a scaling law. Theoretical arguments based upon 

infrared renormalons would favour, instead, a 1/q 2 behaviour. Because of the lack of precise D* production data 
in the intermediate region, it is difficult, at this point, to discriminate between the two possibilities. We point out, 
however, that, if these non-perturbative corrections involve the coefficient functions, they may be present also in light- 
hadron production, where data at intermediate energy are available. It is thus possible that fits to the light-hadron 
fragmentation functions from T(45) up to Z° energies may clarify this issue. 
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